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ABSTRACT 

We present evidence that the very-high-energy (VHE, £ > 100 GeV) gamma-ray emission coincident with 
the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova rem- 
nant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3H-222127, both of which are 
important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected 
emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (a) 

before trials and 7.5 a after trials in a point-source search. The emission is centered at 6*^ 1 6™5 1 ^ + 22°30' 11 " 
(J2000) ±0.03°,,o, ± 0.08°™, with an intrinsic extension of 0.16° ± 0.03°,„„, ± 0.04°„.,. The VHE spectrum 
is well fit by a power law (dN/dE = Nq x (E/TeV)"^) with a photon index of 2.99 ± 0.3Ss,a, ± 0.3,yj and an 
integral flux above 300 GeV of (4.63 ±0.90j,a, ±0.93sv.5) x 10"'^ cm"' s~'. These results are discussed in the 
context of existing models for gamma-ray production in IC 443. 

Subject headings: gamma rays: observations — ISM: individual (IC 443 = VER J0616.9H-2230 = MAGIC 
J0616H-225) 
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1. INTRODUCTION 

IC 443 (G189.1 H-3.0) is one of the most thoroughly stud- 
ied supernova remnants (SNRs) and remains one of the clear- 
est examples of an SNR interacting with molecular clouds. 
IC 443 has a double- shell structure in the optical and radio 
(iBraun & Stroml[r986h . The northeast shell's interaction with 
the H II regi on S249 places IC 443 at a distance of 1 .5 kpc 
( lFesenlll984 b. A second shell is expanding into a less-dense 
region to the southwest. The age of IC 443 remains uncer- 
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tain, with various estimates placing it in the range ~3-30 kyr 
(iPetre et al.ll988HTroja et aLll2008[lLee et alJ2008l:IChevaUeii 

Acc ording to ICornett et alj (Il977h and iDickman et alj 
(11992'), a molecular cloud encircles the remnant. OH 
maser emission from the central and southeast regions in- 
dicates interaction with the foreground p ortion of the cloud 
jClaussen etalJI 19971; iHewitt et al.ll2006i). Th e total mass of 
the clo ud is ^lO" M ^ T ([Torres etal.ll2003h . IDickman etal] 
( 119921) and iLee et alj (|2008|) estimate that -500-2000 Mq 
of cloud material have been directly perturbed by the SNR 
shock. 

IC 443's X-ray emission is primarily th ermal and peaked 
towards the interior of the northeast shell dPetre et 



apeakea 
ill [19881 



iTroia et akl 12006). placing IC 443 in the mixed-morphology 
class of SNRs Observations by Chandra (Olbert et al. 20 01; 
lGaensleretar]20 06) and XMM (iBocchino & Bvkovl200lb of 
the southern edge of the shell have resolved a pulsar wind neb- 
ula (PWN), CXOU J061705.3H-222127, that may be the com- 
pact remnant of IC 443's progenitor. Although the direction 
of motion of the neutron star, inferred from the morphology 
of the X-ray PWN, does not point back to the center of the 
shell, asymmetric expansion of the shell and distortion of the 
PWN tail by turbulence in the ambient medium could be ex - 
aggerating the apparent misalignment (Gaensler et al.ll2006l) . 
No pulsations have been detected at any wavelength. 

The centroid of the unidentified EGRET source 3EG 
J0617H-2238 is located near the center of the IC 443 shell, 
but the limited angular resolution of EGRET makes it im- 
possible to resolve the —45' extension of the remnant com- 
plex. The EGRET source has a spectral index of -2.01 ± 
0.06 over the band 100 MeV - 30 GeV a nd an integral 
flux a bove 100 MeV of 5 x 10"^ cm"^ s"' dHartman et all 
[T999h . AGILE (AGL J0617H-2236) and the Fermi Gamma- 
ray Space Telescope (OFGL J0617.4H-2234) have recently re- 
ported gamma-ray detections with positio n and flux co nsis- 
tent with EGRET (lKttorietal.ll2009l:lAbdoetani2009l ). IC 
443 was detected in very- high-ener gy (VHE, £ > 100 GeV) 
gamma rays by MAGI C d Albert et al. 2007) and VERITAS 
dHumenskv et al.l 120071) . MAGIC J0616H-225 is a point-like 
source coincident with the densest part of the molecular 
cloud and the OH maser emission observed therein. It has 
a power-law spectrum described by (1 .0 ± 0.2^,0, ± 0.35jyj) x 
IQ-^HE/OA TeVr^ 'iO-^-'iO-Z™ TeV"' cm'^ s"' in the range 
0.1-1.6 TeV. In this paper, deep VHE gamma-ray observa- 
tions of IC 443 with VERITAS are reported that reveal ex- 
tended emission coincident with the molecular cloud struc- 
ture. 

2. OBSERVATIONS AND ANALYSIS 



VERITAS dHolder et al.ll2006h consists of four 12-m imag- 
ing atmospheric Cherenkov telescopes located at an altitude 
of 1268 m a.s.l. at the Fred Lawrence Whipple Observatory 
in southern Arizona, USA (31° 40' 30" N, 110° 57' 07" W). 
Each telescope is equipped with a 499-pixel camera of 3.5° 
field of view. The array, completed in the spring of 2007, is 
sensitive to a point source of 1 % of the steady Crab Nebula 
flux above 300 GeV at 5 cr in less than 50 hours at 20° zenith 
angle. 

VERITAS observed IC 443 during two epochs in 2007. 
Data set 1 was acquired with three telescopes during the 
commissioning phase in February and March of 2007. The 

PWN location of 6*^ 1 7™5 .3^ + 22°2 1 '27" (J2000) was tracked 
in wobble mode, in which the source is displaced by 0.5° 



from the center of the field of view dFominI 11994 . Ap- 
proximately equal amounts of data were acquired in each of 
four offset directions. Additional wobble observations with 
four telescopes were taken during October and November of 
2007, using the centroid location determined from the earlier 

data, 6*^16^52.8^ + 22°33'0" (J2000). Application of stan- 
dard quality-selection cuts on weather conditions and hard- 
ware and rate stability resulted in 16.8 hours live time in data 
set 1 and 20.9 hours in data set 2, with mean zenith angles 
of 20° and 17°, respectively. Data set 2 is used for the spec- 
tral analysis. In this data set, 0.8 hours were taken with a 
non-standard array trigger configuration, and these runs are 
excluded from the spectral analysis. 

The da ta are analyzed fol lowi ng standard proced ures de- 
scribed in|Coganiit3] (I1QQ3) andOaniir^^ Show- 
ers are reconstructed for events in which at least two telescope 
images passed pre-selection criteria: number of pixels in the 
camera image > 4, number of photoelectrons in the image 
> 75, and distance of the image centroid from the camera 
center < 1.43°. These criteria impose an energy threshold^^ 
of 300 GeV. 

The event-selection cuts are optimized for weak (—1-5% 
of the Crab Nebula flux) point sources using data taken on 
the Crab Nebula. Cosmic-ray rejection is performed us- 
ing the mean-scaled width (selecting events within the range 
0.05-1.24) and mean-scaled length (selecting events within 
the ran ge 0.05-1.40) parameters, as d escribed in [ Konopelkd 
(119951) . The ring-background model (lAharonian et al.ii2005l) 
is used to st udy the morphology and the reflected-region 
model (lAharonian et al.ll2001l) is used to determine the spec- 
trum. Regions of radius 0.4° around the source location and 
0.3° around two bright stars are excluded from the back- 
ground estimation. All results presented here have been veri- 
fied by an independent calibration and analysis chain; in par- 
ticular, the morphol ogy was cross -checked using a template 
background model dRowelll [20031) because it is sensitive to 
a different combination of systematic effects than the ring- 
background model; this is especially relevant for fields that 
contain bright stars. 

One systematic issue particular to this data set is the pres- 
ence of two bright stars in the IC 443 field: Eta Gem (V band 
magnitude 3.31), located 0.53° from the PWN; and Mu Gem 
(V band magnitude 2.87), located 1.37° from the PWN. The 
high flux of optical photons from the stars requires that sev- 
eral photomultiplier tubes (PMTs) in each camera be turned 
off during the observations and produces higher noise levels 
in the signals of nearby PMTs. These effects reduce the ex- 
posure in the vicinity of the stars and degrade the angular res- 
olution. A high telescope-multiplicity requirement, described 
below, was chosen for studying morphology in order to make 
the direction reconstruction more robust against these effects. 
In order to verify that the presence of a bright star in the field 
does not produce a false excess, a 4.1-hr exposure was taken 
on a sky field that contained a magnitude-2.7 star (36 Eps 
Boo) but lacked any likely gamma-ray sources. No excess 
is observed at the location equivalent to IC 443 in the 36 Eps 
Boo field. 

3. RESULTS 

Here an updated, standard point-source analysis of the 
combined data set (1 & 2) is presented, in which a maxi- 
mum significance of 8.3 standard deviations (cr) before tri- 

■^^ Peak of the differential counting rate for a Crab Nebula-like spectrum. 
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als (7.5 cr post trials, accounting for a blind search over 
the region enclosed by the shell of IC 443) is found at 

06'^16™49^+22°28'30" (J2000) and a significance of 6.8 a 
at the location of MAGIC J06 16+225 . The significan ce is 
calculated according to equation 17 of iLi & Mai (Il983h . Ta- 
ble [T] summarizes the results of this analysis at the location 
of maximum significance, listing the counts falling within the 
point-source integration radius of 0.112° (on), the counts in- 
tegrated in a background ring spanning radii 0.6-0.8° {of f), 
the ratio of on exposure to off exposure (a), and the resulting 
number of excess counts and significance. 

3.1. Morphology 

Figure[T]shows the significance map for the IC 443 field. To 
study the source morphology, an integration radius of 0.112° 
is used and the best-reconstructed gamma rays are selected 
by requiring that an event must have all 3 images and all 
4 images surviving the pre-selection cuts in data sets 1 and 
2, respectively. This requirement increases the analysis en- 
ergy threshold by ^15% and reduces the excess by ~40%. 
The centroid and intrinsic extension of the excess are charac- 
terized by fitting an azimuthally symmetric two-dimensional 
Gaussian, convolved with the PSF^*"* of the instrument, to an 
acceptance-corrected uncorrected excess map with a bin size 
of 0.05°. The PSF has a 68% containment radius of 0.11°. 
The centi-oid is located at 06'^16™51^+22°30'1 1" (J2000) ± 
0.03°j,fl, ±0.08°jv.5, consistent with the MAGIC position. The 
extension derived in this fashion is 0. 16° ± 0.03°^,^, ±0.04°,vv. 
The difference in exten sion between this w ork and the point- 
like source detected by lAlbert et"an (l2007h can be explained 
by the difference in angular resolution and sensitivity between 
VERITAS and MAGIC, the latter of which has an angulai- 
resolution of ~0.15° for 68% containment and a sensitivity 
above 200 GeV of 2% of the Crab Nebula flux in 50 hours 
(lAlbert et alJlIOOSh . 

3.2. Spectrum 

The threshold for the spectral analysis is 300 GeV; the 
energy resolution is less than ~20% at 1 TeV. The pho- 
ton spectrum, integrated within a radius of 0.235°, is shown 
in Figure |2] The photon spectrum is well fit (x^/ndf = 
3 . 1 /3) by a power law dN/dE =Nox{E /TeV)-^ in the range 
0.317-2.0 TeV, with a normalization of (8.38 ± 2.10,,fl, ± 
2.50,ys) X 10"'^ TeV"' cm'^ s"' and an index of 2.99 ± 
0.38,,fl, ±0.30,vs. The integral flux above 300 GeV is (4.63 ± 
0.90,,fl,±0.93,vj) X 10-12 s"' (3.2% of the Crab Nebula 
flux), consistent within errors with the spectrum reported by 
MAGIC (.Albert et al.. 2007). 

4. DISCUSSION AND CONCLUSIONS 

Figure [3] shows the inner 0.8° of the excess map and places 
the VHE emission in a multi-wavelength context. The TeV 
centroid is ~ 0. 15° from the position of the PWN and ~ 0.03° 
from the nearby, bright m a ser em ission coincident with clump 
G of iHuang & ThaddeusI ( 119861) . Two additional weaker re- 
gions of maser emission lie along the souther n rim of the 
remnant, east of the PWN dHewitt et al.ll2006h . The VHE 
emission overlaps the foreground molecular cloud. While the 

■^"^ The PSF is characterized as a sum of two, two-dimensional Gaussians 
describing a narrow core and a broader tail. It is determined from data taken 
on th e Crab Nebula, which is a point source at these energies lAlbert et alj 
l200l . 



Fermi source OFGL J06 17. 4-1-2234 is displaced from the cen- 
troid of the VHE emission by ^0.15°, this is consistent at 
the 95% level with the combined errors quoted by Fermi and 
VERITAS. 

The VHE gamma-ray emission observed by VERITAS and 
MAGIC is offset from the location of the PWN by -10-20 
arcmin, co mmensurate with other of fset PWNe such as HESS 
J1825-137 dAharonian et al.ll2006al) . The emission is con- 
sistent with a scenario in which the VHE emission arises 
from inverse Compton scattering off electrons accelerated 
ear ly in the PWN's li fe. Such a scenario was presented 
by iBartko & Bednare^ ((2008). Under the assumption that 
the VHE emission is associated with IC 443 at a distance 
of 1.5 kpc, the luminosity in the energy band 0.3-2.0 TeV 
is 4 X 10^2 erg s"'. The spin-down lum inosity of the pul- 
sar has been estimated as '^lO^^^ erg s"' (lOlbe rt et al.ll200ll: 
Bocchino & Bykov 2001) and 5x10 erg s"' (Gaensle r et alj 
l2006t) . If the PWN association is correct, the VHE luminosity 
in the 0.3-2.0 TeV band is less than —0.04% of the spin-down 
luminosity, well within the —0.01-10% range of observed 
VHE efficiencies for other PWNe (Gallant 2007). Likewise, 
synchrotron cooling is unlikely to rule out an offset PWN 
scenario. The Ufetime r of a synchrotron-emitting electron 
wifli energy £ is t(E ) - 1.3 x lO'^(B/^iG)-^{E/l GeV)"' kyr 
dGaisser et al.l Il998l) . Assuming the TeV photons are pro- 
duced by — 20-TeV electrons in a magnetic field of —5 fiG 
(comparable to typical interstellar fields), synchrotron cooling 
is only just beginning to become important even if the age of 
IC 443 is as high as —3 kyr. Note that in the nebu la-powered 
scenario presented by 'B artko & Bednare^ (l2008h . the GeV 
emission is assumed to be unresolved pulsar emission. The 
displacement of OFGL J0617.4H-2234 from the PWN location 
argues against that scenario. 

Alternatively, Figure[3]can be interpreted within a scenario 
of hadronic cosmic-ray acceleration and subsequent interac- 
tion with the molecular cloud, which would provide a high 
density of target material for the production of VHE gamma 
rays. The correlation of the VHE emission with the molec- 
ular cloud is natural within this scenario. The low velocity 
of the SNR shock implied by the presence of maser emis- 
sion (Hewitt et al. 2006) indicates that the shock is in its ra- 
diative phase in this region and cosmic-ray acceleration is 
most likely now inefficient. Thus, the cosmic rays acceler- 
ated during the shock's earlier propagation have diffused into 
the remnant and the upstream region, including the molec- 
ular cloud. The steep VHE spectrum can be explained ei- 
ther as a low maximum energy to which particles were ac- 
celerated prior to the shock hitting the cloud, or an energy- 
dependent rate of diffusion of cosmic rays out of th e cloud 
dAharonian & AtoyanI [T996h . IZhang&FangI d2008h model 
the remnant as evolving partially within and partially outside 
the cloud, and find that a hadr onic interpret ation can explain 
the VHE e mission. However. [Torres et alj f2008) point out 
that in flie IZhang & Fanj d2Q08|) model, flie GeV and TeV 
emission ought to be spatially coincid ent; it is not yet clear 
whether they are. [Torres etakl d2008l) present an alternative 
picture in which cosmic rays escape from the SNR shock and 
diffuse toward the cloud, similar to what may be occurring 
in the Galactic Ridge ( Aharonian et al. 2006b) and the region 
surrounding the SNR W 28 (A haronian et al.a2008i) . This pic- 
ture can naturally accommodate spatially separated GeV and 
TeV emission b y positing clouds at diffe rent distances from 
the SNR shock. iRodriguez Marrero eFal. (2009) extend this 
work to make predictions for the Fermi Gamma-ray Space 
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Telescope (Fermi) and find that Fermi should see a movement 
of the centroid towards the TeV position as the photon en- 
ergy increases, along with a gradual spectral softening. The 
combination of GeV and TeV observations may also provide 
additional constraints on the diffusion coefficient and the dis- 
tance between the SNR and the molecular cloud. 

Further GeV and TeV gamma-ray observations will help 
to clarify the nature of the particle acceleration associated 
with IC 443. The ^0.15° separation between the centroids 
of OFGL J0617.4H-2234 and VER J0616.9H-2230, though not 
statistically significant, may hint at an energy-dependent mor- 
phology. Such a morphology will be explored with future ob- 
servations by VERITAS, and has the potential to discriminate 
between the scenarios outlined above. A break in the energy 
spectrum between the GeV and TeV bands is required whether 
or not OFGL J0617.4H-2234 and VER J0616.9H-2230 are di- 
rectly related. Future observations by VERITAS will improve 
the statistics above 100 GeV, anchoring the measurement of 
the VHE cut-off to the gamma-ray emission. The combined 
spectrum from 100 MeV to several TeV will provide power- 
ful co nstraints on the emission mechanism(s) (G aisser et alj 
fT998h . 

In summary, VERITAS observations of the composite su- 
pernova remnant IC 443 during 2007 have yielded a detection 



of an extended source with an index of 2 . 99 ± . 3 8 ,,0, ± . 30jvs 
and aflux of (4.63±0.90j,n,±0.93„.,) x 10"'- cm'^ s"' above 
300 GeV. The location and flux are consistent with MAGIC 
J0616H-225. This deep VERITAS observation reveals that 
the VHE emission is extended (0.16° ± 0.03°,,<„ ± 0.04°,„), 
with its brightest region coincident with the dense cloud ma- 
terial and maser emission. The emission is offset from the 
nearby PWN CXOU J061705.3H-222127, a viable source for 
the VHE emission. If further VHE observations reveal an 
energy-dependent morphology, it may become possible to dis- 
tinguish cleanly between scenarios related to the PWN and to 
hadronic cosmic rays interacting with the molecular cloud. 
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Fig. 1. — Significance map for the IC 443 field. The black circle indicates 
the one-sigma angular extension (see text). The black cross-hair indicates 
the centroid position and its uncertainty (statistical and systematic added 
in quadrature), and the white cross-hai r likewise ind icates the position and 
uncertainty of MAGIC J0616+225 ( Albert e t alJ20Q7D . Red contours: optical 
intensity ( McLean et al. 2000). White star: location of the star Eta Gem. The 
white circle indicates the PSF of the VERITAS array. 



TABLE 1 

Analysis results for VERITAS observations of IC 443. 





on 


off 


a 


excess 


significance 
(cr) 


data set 1 


393 


5152 


0.0610 


78.7 


4.1 


data set 2 


609 


7331 


0.0601 


168.1 


7.3 


combined 


1002 


12483 


0.0604 


247.5 


8.3 
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Fig. 2.— Spectram of IC 443, scaled by E'. Red points are the VERITAS 
spectrum (red tick marks along top indicate bin edges). The red line is a 
power-law fit (see text), with residuals to the fit plotted in the lower box. The 
MAGIC spectrum is indicated by the gray band and extends down to 90 GeV. 
VERITAS eiTor bars and MAGIC error band reflect statistical en'ors only. The 
Crab spectrum (V. Acciari et al. 2009, in preparation) is shown as a dashed 
line for comparison. 
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Fig. 3. — Inner 0.8° of the acceptance-corrected excess map for the IC 
443 field. Markers and contours follow Figure [T ] with sev eral additions. 
Thick black contours: CO survey ( Huans & Thaddeus 1986); black star: 
PWN CXOU J061705.3+222127 (Olbert et al. 2001); open blue circle: 95% 
confidence radius of OFGL J0617.4+2234 tAbdo et al..2009,) ; and filled black 
triangles: locations of OH maser emission (Claussen et al. 1997; J. W. Hewitt, 
private communication). 



